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THRESHOLD  CRITERIA  FOR  A  SPACE-SIMULATION  BEAM-PLASMA-DISCHARGE 

I.  INTRODUCTION 

A  cold  electron  beam,  propagating  through  a  weakly 
ionized  plasma  will,  under  proper  conditions,  produce  a 
modified  beam-plasma  state  known  as  the  Beam-Plasma-Discharge 
(BPD).  This  discharge  state  has  received  considerable  attention 
in  recent  years  as  a  result  of  increased  interest  in  mechanisms 
for  vehicle  neutralization  during  spaceborne  accelerator 
experiments  (Bernstein,  et  al . ,  1980;  Cambou,  et  al  .  ,  1978), 
enhanced  beam-plasma  ionization  processes  (Bernstein,  et 
al.,  1978),  and  in  general  single-particle  or  collective 
phenomena  initiated  by  beams  injected  into  neutral  gas  and 
charged-particle  environments  (Hess  et  al.,  1971;  Winckler, 
et  al.,  1975;  Hendrickson  and  Winckler  1976;  Cambou,  et  al., 
1975;  Munson  and  Kellogg  1978a;  Szuszczewicz ,  et  al.,  1979; 

Jost  et  al.,  1980;  Winckler  1980). 

The  BPD  appears  at  a  critical  energetic-electron-beam 

c 

current  I_,  with  the  transition  from  single-particle  behavior 

D 

(I  <  1^,  pre-BPD)  to  collective  processes  (I  >  Ic,  solid- 
B  B  B  B 

BPD)  described  as  follows  for  conditions  in  which  the  plasma 
is  created  by  the  beam  itself: 

(i)  As  an  electron  beam  linearly  interacts  with  a 
neutral  gas,  it  c o 1 1 i s  i  on a 1 1 y  produces  a  plasma  with  a 
density  that  varies  directly  with  the  magnitude  of  the  beam 
current  for  a  fixed  beam  energy. 

(ii)  As  the  beam  current  is  increased  to  a  critical 
c 

value,  1^,  a  two-stream  instability  sets  in  and  the  electric 
fields  of  the  excited  waves  heat  the  electrons  to  energies 
comparable  to  the  ionization  energy.  The  "heated"  electrons 
create  an  enhanced  ionization  process  which  results  in  an 
avalanche  breakdown,  the  BPD. 

Manuscript  submitted  July  14.  1981. 
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Bernstein  et  al.,  (1979)  have  reported  the  dependence 

c 

of  this  critical  current,  Ifi,  on  various  experimental 

parameters  as 

V  3/2 

— _  (1) 

B  b°'7pl 

where  V_  ,  B,  P  and  L  are  the  beam  energy  (voltage),  the 

D 

superimposed  magnetic  field,  the  ambient  neutral  pressure, 
and  the  beam  length  (gun-to-collect or  distance),  respectively. 
While  the  (V_,B,P,L)  relationship  was  established 

D  D  D 

among  the  controlling  system  parameters,  a  clear  dependence 

on  plasma  density  was  expected,  with  early  thoughts  suggesting 

that  a)  -  a)  satisfied  ignition  threshold  criteria  (Bernstein, 
P  c 

et  al.,  1979;  Getty  and  Smullin,  1963).  We  have  since  had 
the  opportunity  to  test  these  ideas  under  various  beam- 
plasma  conditions.  Our  results  include: 

(i)  Time-delay  data  and  associated  analyses  indicating 
that  BPD  ignition  does  indeed  occur  at  a  critical  density, 

(ii)  Direct  measurements  of  plasma  density  near  ignition 

threshold  with  determinations  of  0)  / to  ,  and 

P  c 

(iii)  A  theoretical  analysis  which  predicts  the  critical 
density  criterion. 


In  subsequent  sections,  we  present  the  experimental  and 
theoretical  details  which  establish  the  density-dependent 
threshold  conditions  for  BPD  at  3.9  £  to  /to  £  5.8. 

K  C 

II.  EXPERIMENT  CONFIGURATION  AND  RESULTS 


The  experiment  was  conducted  in  a  large  vacuum  chamber 
facility  at  the  NASA  Johnson  Space  Flight  Center  with  an 
experimental  configuration  similar  to  those  employed  in 
earlier  investigations  (Bernstein  et  al.,  1979).  The 
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configuration  is  illustrated  in  Figure  1  with  specific 
emphasis  on  the  pulsed  plasma  probe  measurement  and  the 
associated  procedure  for  oo  /u)^  determinations.  Not  shown 
are  the  previously  described  diagnostics  (Bernstein  et  a  1  .  , 

1979;  Jost  et  aL ,  1980)  including  the  39148  scanning  photometer, 

the  segmented  current  collector,  the  energetic  electron 
electrostatic  analyzer,  and  the  remote  wave  detection  antenna 
and  spectrum  analyzer  system.  The  beam  was  generated  by  a 
tungsten-cathode  Pierce-type  diode  gun,  mounted  on  a  position- 
controlled  cart.  In  all  cases  the  beam  was  injected  parallel 
to  the  magnetic  field  B  and  terminated  on  a  3  x  3  m  target 
suspended  about  20  m  above  the  gun  aperture.  A  combination 
of  coil  current  and  the  Earth's  magnetic  field  established 
the  B-field  at  selected  levels  up  to  1.5  gauss. 

In  most  cases  the  beam  was  injected  into  a  neutral  gas 
with  no  pre-beam  plasma;  however  the  experimental  investigation 
included  two  cases  in  which  the  chamber  was  filled  with  a 
plasma  created  by  a  Kauffman -type  argon  ion  thruster.  In 
these  cases  the  pre-beam  plasma  density  was  lower  than  the 
critical  density  at  BPD  ignition. 

The  investigation  was  conducted  in  two  stages.  First, 
pulsed-gun  experiments  were  carried  out  to  test  the  concept 
of  a  d e n s  i  t y -d e p e nd en t  threshold  condition  against  the 
original  ideas  of  Getty  and  Smullin  (1963).  After  the 
concepts  check,  a  series  of  direct  density  measurements  were 
conducted  to  quantitatively  establish  the  density  criteria. 

We  describe  the  results  in  the  order  in  which  the  experiments 
were  conducted. 


c 

Time-Delay  Arguments  for  N^ 

_  c 

In  pulsed-beam  experiments  with  I  >  >  I,  Getty  and 

D  D 

Smullin  (1963)  identified  three  sequential  phases  in  the 
temporal  evolution  into  the  BPD:  1)  a  quiescent  stable 
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Fig.  1  —  Experiment  configuration 


period,  2)  an  intermediate  period  during  which  the  gross 

features  of  the  beam  remained  unchanged  but  waves  with  f  fv 

f  were  present,  and  3)  BPD  ignition.  They  suggested  that 
c  e 

phases  1  and  2  corresponded  to  the  time  required  for  buildup 

of  the  ambient  density  to  a  critical  value. ..a  value  at 

which  the  BPD  ignited.  To  test  this  concept  in  the  space- 

simulation  beam-plasma  efforts  at  the  JSC  facility  a  series 

of  pulsed  beam  measurements  were  performed  to  study  the 

temporal  evolution  of  the  BPD  in  the  configuration  shown  in 

Fig.  1.  All  diagnostics  sensitive  to  the  transition  to  BPD 

demonstrated  identical  time  delays  between  beam  current 

_ b 

initiation  and  BPD  ignition.  At  low  pressure  (<  4(10  ) 

torr)  the  Getty  and  Smulli n  phase  2  (characterized  by 

f  %  f  waves)  was  clearly  evident.  Figure  2  shows  the 
ce 

measured  delay  times  as  a  function  of  beam  current  with 

values  ranging  from  0.1-20  msec.  The  results  in  Figure  2 

can  be  understood  in  terms  of  a  simple  time-dependent  model. 

If  it  is  assumed  that  the  plasma  loss  rate  is  proportional 

to  its  density,  then  the  temporal  buildup  of  plasma  by 

electron-neutral  collisions  is  given  by 

N  =  —  *  (2) 
e  ^  (1  -  exp  [  -  —  t  ] ) . 


t  ]). 


Equation  (2)  follows  directly  from  Getty  and  Smullin  (1963) 
as  the  solution  to 


=  P  -  ~  N 
L  e 


where  P  =  In  ON  /eA  is  the  electron-ion  pair  production 
Bo 

rate,  A  /  L  is  the  loss  rate,  O  the  electron  impact  ionization 

cross  section,  N  the  neutral  density  and  A  the  beam  cross 
o 

section.  At  short  times,  the  density  will  increase  linearly 
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BEAM  CURRENT 


Legend-'  a)  Measured  Delays  • 
b)  Theoretical  Fit  - 


Conditions1  E  =  1.8  keV,  B  =  1.56  gauss 
P-  6(10”*)  torr 

WToff*  80  ms/250  ms 


BPD  Threshold  Current 


0  20  40  60  80  100  120  140  160  180 

BEAM  CURRENT  [ma] 

Fig,  2  —  Time  delay  in  the  ignition  of  the  beam-plasma-discharge 
with  threshold  current  at  20  ma 
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with  time  (N^  =  Pt)  and  at  longer  times  the  density  will 
approach  the  time-independent  value  =  PL/A.  For  1^  ^ 

Q 

I  relatively  long  delay  times  should  be  observed  with  a 
B 

rapid  decrease  in  t  as  I  is  increased.  This  is  exactly 

a  d 

the  behavior  demonstrated  in  Figure  2,  where  Equation  (2) 
has  been  normalized  to  the  data.  The  good  agreement  between 
the  simple  model  calculations  and  the  experimental  results 
supports  the  concept  of  a  critical  density  threshold  for  BPD 
ignition,  but  the  normalization  procedure  provides  little 
measure  for  an  absolute  value  of  N^.  Specific  determinations 
of  the  critical  density  threshold  criteria  are  detailed  in 
the  next  two  sections. 


Direct  Measurements  of  Critical  Density 

To  quantitatively  establish  the  density-dependent 

threshold  criterion  for  the  ignition  of  a  space-simulation 

BPD,  emphasis  was  placed  on  the  direct  measurement  of  plasma 

density  profiles  over  a  range  of  beam-plasma  conditions 

covering  beam  currents  from  8  to  85  ma,  beam  energies  from 

0.8  to  2.0  keV  and  magnetic  fields  at  0.9  and  1.5  gauss. 

The  procedure  involved  experimental  determinations  of  radial 

profiles  of  electron  density  for  each  of  the  selected  conditions 

extending  from  a  low-pre-BPD  state  to  a  strong  BPD  condition. 

The  experimental  configuration,  illustrated  in  Figure  1, 

involved  a  pair  of  pulsed-plasma-probes  mounted  on  a  radial 

traversal  mechanism  positioned  at  approximately  8  m  above 

the  injection  point  of  the  beam.  Each  of  the  probes  provided 

simultaneous  measurements  of  eLectron  density  N^,  temperature 

,  plasma  potential  V m  ,  and  density  fluctuation  power 

spectra  6N  (+P  (k))  with  capabilities  for  the  associated 

e  n 

diagnostics  in  a  dynamic  plasma  environment  and  under  test 
conditions  which  could  contaminate  electrode  surfaces  (Holmes 
and  Szuszczewicz,  1975,  1981;  Szuszczewicz  and  Holmes  1975,1976). 


The  plasma  density  measurements  were  made  for  seven 

different  conditions,  each  identified  by  pre-selected  values 

for  VD,  B,  P  and  the  existence  or  non-existence  of  a  pre- 
B 

beam  plasma.  For  each  condition  a  steady  state  value  for  I 

B 

was  set,  the  traversal  mechanism  exercised,  and  an  electron 
density  profile  was  recorded.  A  sample  profile  collected 
under  pre-BPD  conditions,  is  presented  in  Figure  3.  The 
abscissa  is  time  relative  to  the  start  of  the  radial  traversal 
and  the  ordinate  is  relative  electron  density  as  determined 
by  baseline  electron-saturation  currents  collected  by  the  E - 
probe.  (The  second  in  the  two-probe  configuration  was 
defined  as  the  1-probe  because  the  associated  baseline 
currents  were  collected  in  the  ion-saturation  portion  of  the 
probe* s  current-voltage  characteristic  (Holmes  and  Szuszczewicz, 
1973,  1981).)  At  the  start  of  each  traversal  the  probe  was 

at  its  outermost  position  relative  to  the  center  of  the 
chamber.  As  time  increased  the  probe  was  moved  into  and 
through  the  beam;  at  maximum  penetration  toward  the  chamber 
center,  the  traversal  system  was  reversed,  allowing  a  second 
measurement  of  the  density  profile  as  the  probe  moved  back 
to  its  original  outermost  position.  With  this  procedure  the 
probe's  maximum  penetration  toward  the  chamber  center  is 
identified  by  the  symmetry  point  in  the  "double"  profile. 

The  symmetry  in  the  "double"  profile  provided  confidence 
that  beam-plasma  conditions  were  unchanged  during  the  measure- 
men  t . 

Absolute  electron  densities  were  determined  by  standard 

analysis  procedures  summarized  graphically  in  Figure  4. 

3 

The  technique  provides  a  high-frequency  (10  Hz  in  this 
experiment)  determination  of  relative  electron  density 
through  the  direct  measurement  of  baseline  e 1 e c t r on - s a t u r a t i on- 
currents.  Simultaneously,  the  technique  generates  a  "conventional" 
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t  J, 


tana* 


TIME 


Fig.  3  —  Radial  profile  of  relative  electron  density  under  pre-BPD  conditions.  Run 
#57,  (IB,  VB,  B)  =  (7  ma,  1.3  keV,  0.9g).  The  figure  shows  two  cuts  through  the 
beam-plasma  profile,  as  time  increases  from  left-to-right  the  plasma  density  probe 
moves  into  and  through  the  beam  center,  then  reverses  and  passes  through  the  beam  a 
second  time.  The  symmetry  verifies  that  beam-plasma  conditions  were  stable  during 
the  execution  of  the  radial  traversal. 
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Fig.  4  —  Sample  of  raw  probe  data  (3A)  showing  the  effects  of  density 
fluctuations  (baseline  electron-saturation-currents)  on  the  probe’s  current- 
voltage  characteristic  (sweep  currents).  3B  shows  the  “corrected  ”  charac¬ 
teristic  that  results  when  density  fluctuations  have  been  unfolded.  The 
electron  saturation  portion  of  the  characteristic  is  then  amenable  to  con¬ 
ventional  analyses  for  determination  of  Ne  (see  text). 
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Langmuir  probe  characteristic.  The  relative  density  fluctuations 
(as  indicated  by  the  variations  in  the  baseline  current)  are 
then  unfolded  from  the  raw,  uncorrected  probe  characteristic 
(Fig.  4A)  yielding  a  smooth,  corrected  curve  (Fig.  4B)  to 
which  conventional  N^-analysis  procedures  (Chen,  1965; 

Szuszczewicz  and  Holmes,  1977)  are  applied.  This  procedure 
was  utilized  for  all  beam-plasma  conditions  included  in  this 
inves  t igation . 

Complete  profile  information  and  associated  plasma  wave 
signatures  are  presented  for  two  independent  conditions  in 
Figures  5  and  6.  (Note:  Because  of  changes  in  plasma 
potential,  the  relative  density  profiles  do  not  maintain  the 
same  scaling  to  absolute  values  from  run-to-run.)  For 
values  of  (Vb,B)  =  (1.3  keV,  0.9G)  and  (2.0  keV,  1.5G)  in 

D 

Figures  5  and  6,  respectively,  the  beam  current  I  was 

B 

stepped  through  a  sequence  allowing  for  complete  coverage  of 

conditions  which  encompassed  pre-,  threshold-,  and  solid- 
£ 

BPD.  Because  1  D  =  I D  ( V  D  ,  B  ,  P  ,  L  )  is  known  to  exhibit  hysteresis 

D  D  D 

as  lD  is  varied  about  the  critical  value  I  ^  ,  the  following 

D  D 

procedure  was  utilized  in  the  step-wise  selection  of  beam 
current  levels: 

(i)  With  V  fixed,  L  was  steadily  increased  in  a 

D  O 

"search"  mode  to  determine  an  approximate  value  for  1  . 

B 

c 

(ii)  I  was  reset  to  zero,  then  increased  slowly  to  I 
d  B 

with  careful  observation  of  the  RF  spectrum.  Threshold  was 
defined  as  that  level  at  which  the  RF  spectrum  was  in  transition 
from  its  pre-  BPD  signature  (flat  spectrum  at  high  frequencies, 
f  f  )  to  its  solid-BPD  characteristic  (intense  features  at 

>  P 

f  ^  f  ).  At  this  transition  level  (defined  here  as  threshold) 

P 

the  beam-plasma  system  oscillated  between  its  two  states  at 
a  repetition  frequency  generally  observed  to  be  no  greater 
than  several  Hertz.  If  the  RF  spectrum  "locked-up"  in  a 
characteristically  stable  RF  BPD  signature,  step  "(ii)" 
was  restarted  and  the  previous  1^  level  more  carefully  approached. 
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RELATIVE  N„  RELATIVE  N 


V\£.  5  Sequence  of  relative  plasma  density  profiles  ( h’114.  IA)  and 
associated  plasma  wave  signatures  (Fi^.  *1B|  for  increasing  values  of  beam 
currenl  1^  at  a  fixed  condition  |V^,B)  (1.3  keV,  O.tyj)  encompassing 

runs  5ft  through  83  (pre-BI*l>  through  «olid-BI*D). 


(b) 

Fig.  B  —  Sequence  of  relative  plasma  density  profiles  (Fig.  5A)  and 
associated  plasma  wave  signatures  (Fig.  5B)  for  increasing  values  of  beam 
current  l^  at  a  fixed  condition  (V^B)  ~  (2.0  keV,  1.5g)  encompassing  runs 
79  through  83  (pre-BPI)  through  solid-BPD). 


(It  is  important  to  note  that  the  time-dependent  photographs 

of  the  threshold  RF  spectrum,  e.g.,  runs  57  and  81,  erroneously 

indicate  a  solid-BPD  spectrum  and  do  not  reflect  the  oscillations 

between  the  two  beam-plasma  states.)  The  conditions  at 

threshold  and  under  BPD  are  summarized  in  Table  1  where  the 

peak  density  NmaX,  associated  plasma  frequency  2TTU)maX, 

max  ^ 

and  plasma-to-cyclotron  frequency  ratio  a)  /oj  are  also 

P  C  +1  3 

listed.  The  results  can  be  summarized  by  u>  /u>  *  5.8  *  as 

pc  -1.9 

the  density-dependent  threshold  condition  for  the  BPD.  As 

we  will  show  in  the  next  section,  this  result  is  consistent 

with  a  threshold  model  which  assumes  that  BPD  is  triggered 

by  the  onset  of  a  beam  plasma  instability. 

It  is  proper  to  introduce  one  qualification  regarding 

the  experimentally  determined  values  for  w^/u)  *  By  definition, 

P  c 

they  represent  averages  between  two  beam-plasma  states 

(pre-and  solid  BPD).  The  consequence  of  this  in  the  interpretation 

of  the  experimental  values  is  merely  one  of  definition.  For 

example,  if  threshold  were  redefined  as  the  maximum  density 

under  stable  pre-BPD  conditions  (i.e.,  no  oscillations  into 

the  solid-BPD  state),  then  (jo^/o)  =  5.8  would  represent  an 

upper  limit.  Our  experience  suggest  that  the  "ad j ustment M  in 

critical  density  would  not  be  more  than  a  factor  of  two  less 

than  the  measured  values.  In  this  case  we  suggest  that  a 

more  general  statement  of  density-dependent  threshold  take 

the  form 

3.9  'V  u>  /w  <  5.8.  (A) 

P  c 

Before  proceeding  to  the  theoretical  treatment  it  is 
appropriate  to  discuss  several  of  the  plasma  density  and 
wave-related  features  in  the  data  presented  in  Figures  5  and 
6. 


Run  56  (upper  left  panel  Figs.  5A  and  5B)  is  very 
typical  of  a  pre-BPD  environment,  in  that  the  team  center 
(see  e.g..  Fig.  3)  is  well  defined  and  of  relatively  narrow 
cross  section,  with  an  RF  signature  (upper  left  panel  Fig. 
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5B)  that  is  flat  in  character.  At  BPD  threshold  the  beam- 
plasma  system  intermittently  emits  RF  at  and  around  the 
plasma  frequency  and  the  plasma  expands  in  cross  section 
(lower  left  panel  Fig.  5A)  eliminating  the  narrow  beam 
profile  characteristic  of  the  pre-BPD  state  (Fig.  3  and  5A 
upper  left).  We  note  that  the  RF  signatures  are  highly  time 
dependent,  generally  fluctuating  in  intensity  and  in  frequency. 
Similarly,  the  local  plasma  densities  are  highly  time  dependent 
(see  e.g.,  fluctuations  in  Fig.  4)  with  that  dependency 

averaged  over  i-sec  intervals  for  the  relative  density 
profiles  in  Figs.  5  and  6. 

As  I  is  increased  above  the  threshold  level  (Runs  b  b 
ii 

and  59  in  Figs.  5A  and  5B)  the  plasma  density  increases  and 
the  radial  profiles  take-on  variable  configurations.  In 
Run  58,  for  example,  the  beam  center  is  "depleted"  with  an 
increased  density  level  at  its  edges.  This  configuration 
has  been  observed  in  a  number  of  cases  and  is  thought  to 
manifest  electron  heating  in  the  beam  core  with  generally 
higher  diffusion  rates  and  associated  losses.  The  profile 
becomes  more  complex  in  Run  59,  possibly  due  to  higher  order 
diffusion  modes. 

The  results  in  Figure  6  have  been  selected  to  present 
information  at  a  higher  applied  magnetic  field  and  to  illustrate 
that  the  description  of  beam-plasma  profiles  offered  in 
connection  with  Figure  5  is  by  no  means  considered  universally 
applicable  to  the  various  levels  of  the  BPD.  Our  current 
level  of  understanding  does  not  allow  for  a  detailed  cause- 
effect  description  of  each  and  every  profile  but  at  present 
we  would  point  to  Figure  5A  as  the  simplest  configuration 
that  illustrates  the  transition  from  pre-  BPD  to  the  solid 
BPD  state. 

III.  THEORETICAL  CONSIDERATIONS  AND  DISCUSSION  OF  RESULTS 

Before  presenting  the  theoretical  details  of  the  BPD 
threshold  criteria  as  applied  to  the  present  experiment,  it 


is  important  to  discuss  some  aspects  of  our  model  which  play 
a  key  role  in  the  interpretation  of  the  experimental  results. 

For  all  experimental  parameters  the  relevant  mean  free 
paths  (i.e.,  thermalization,  ionization,  etc.)  were  much 
longer  than  the  system  length  L.  This,  as  discussed  in 
detail  by  Papadopoulos  (1981)  for  a  finite  system,  leads  to 
a  requirement  for  axial  confinement  of  both  the  ionizing  and 
the  ambient  electrons.  As  shown  in  the  above  paper,  the 
required  axial  confinement  can  be  achieved,  if  the  energy 
deposited  by  the  beam  electrons  is  absorbed  by  few  ambient 
electrons  (i.e.,  1%)  which  are  accelerated  to  energies 

larger  than  the  ionization  potential.  A  fraction  of  these 
electrons  quickly  escapes  the  system,  thereby  building 
potential  sheaths  at  the  ends.  The  sheaths  accelerate 
ambient  ions  to  escape  at  the  rate  of  the  fast  electrons, 
while  trapping  the  cold  part  of  the  energetic  electron 
population.  A  detailed  description  of  the  process  can  be 
found  in  Papadopoulos  (1981).  For  the  purposes  of  the 
present  paper  we  retain  the  requirement  that  the  collisionless 
mechanism  responsible  for  BPD  should  deposit  most  of  its 
energy  to  suprathermal  tails  with  energy  much  larger  than 
the  ionization  energy  while  leaving  the  majority  of  the 
electrons  cold.  As  discussed  in  many  publications  (Papadopoulo 
and  Coffey  1974;  Papadopoulos  and  Rowland  1978;  Linson  and 
Papadopoulos  1980;  Papadopoulos  1981),  this  can  he  achieved 

If  the  frequence  of  the  excited  waves  is  near  the  electron 

W 

plasma  frequency  (i>  ,  and  the  wave-energy  — — -  is  greater 

2  2  P  N  1 

than  k  where  W  is  the  field  energy  density,  N  the 

ambient  electron  density,  T^  the  electron  temperature,  k  the 

wave-number  of  the  instability  and  A  the  Debye  length.  In 

W  u 

order  for  substantial  to  be  built  in  a  finite  length 

system  the  instability  should  be  almost  absolute.  The 
dispersion  relation  for  a  beam  plasma  system  is  (Linson  and 
Papadopoulos,  1980;  Rowland  et  a  1  .  ,  1981) 


k2(l-  ^£)  +  k2  (1  -  -^~2  ) 

Z  2  i-  co-co 

0)  c 


2  ou)2  2 

kz  -  2  "  k 

U)-kzVb) 


■■  0  (5) 


(u)-k  v  )2-0)2 

ZD  C 


where  k 


with  T  ,  the  roots  of  the  Bessel  function 
s  1 


J  (T  .)  =  0  and  r  the  beam  radius.  In  Eq .  (5),  a  is  the 

1  s  1  o 

b  earn- 1  o- p  1  a  sma  density  ratio,  0)^  and  the  plasma  and 

electron  cyclotron  frequencies,  the  beam  velocity  and 

the  parallel  wave  number.  On  the  basis  of  Eq.  (5)  we  can 

find  (Rowland  et  al . ,  1981)  that  an  almost  absolute  instability 

can  be  excited  near  u)  ,  for 

e 


_£  -  l 

2 


%  p 

Coupling  this  inequality  with  k^  %  ^  and  k  - 

our  required  threshold  condition  as  ^ 

2  2  2 
(2.4)  vu  cosZ0 
„  b 


gives 


The  criterion  (7),  with  0  being  the  beam's  angle  of  injection 
relative  to  the  magnetic  field,  can  be  generalized  to  the  case 
where  the  plasma  radius  R  is  different  from  the  beam  radius  rQ 
(Rowland  et  al.,  1981).  In  this  case  Eq .  (7)  becomes 

2  >  vb  COs2°  (8) 

0)  -  2  — - - 

p  r2  ln(R/r) 

o  o 

The  value  of  r  is  given  by  considering  both  the  Larmor 
o 

radius  as  well  as  the  term  including  space  charge  beam 
expansion  (Linson  and  Papadopoulos ,  1980).  This  gives 

r  -  2.1  B  ,  (In  °  .  sin  20)>S  m  (9) 


where  <  is  gun  perveance  in  micropervs,  rg  the  gun  radius  in 
cm,  IR  and  V  the  beam  current  in  amps  and  energy  in  keV , 
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and  the  magnetic  field  B  is  in  gauss.  From  Eq .  (8)  and  (9) 

%  r  o  'v  1 

and  using  r  2mm  and  — r  ^  *r  as  observed  we  find 

g  R  2 


5x10  1Vr3/4 

D 


co  s  0 


/  ,  i  8*5  .  2  n  v 

(1  +  — ~ —  sin  0 ) 


We  can  see  that  the  u)  /  U)  ratio  at  threshold  (with  k  ^  1 

P  c 

microperv)  is  rather  insensitive  to  the  beam  injection  angle 

for  0  'v  60°.  In  Table  II  we  present  the  values  of  w  /u) 

P  c 

computed  on  the  basis  of  Eq.  (10)  by  taking 
cos  0 

8.5  2  f\j 

(1  +  sin  ®)2  ^  1  ^or  near  parallel  injection.  From 

Table  II  we  find  the  average  computed  value  to  be  u)  /to 

P  c 

4.95.  This  result,  while  subject  to  moderate  uncertainties 
in  r^/R,  is  taken  to  be  in  excellent  agreement  with  the 
experimentally  derived  conditions  (Table  1  and  Eq.  4)  providing 
complementary  arguments  which  confirm  the  original  suggestion 
that  the  ignition  of  the  BP D  was  strongly  coupled  to  a 
density-dependent  threshold  criterion.  We  should  mention 
that  the  above  theory  is  consistent  with  the  observational 
facts  that  the  temperature  of  the  electrons  is  much  lower 
than  the  ionization  potential  (Szuszczewicz  et  al . ,  1979), 

while  energetic  tails  are  observed  within  the  beam-plasma 
core  (Szuszczewicz,  1980  unpublished). 


Run  No. 


TABLE  II:  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 


Vfi(keV) 


(ma) 


u>  /co  (calc.) 
P  c 


U)  /(a) 

P  c 


(observed) 


40 

1.9 

37 

4.2 

6.9 

48 

1.9 

34 

4.4 

6.6 

57 

1.3 

18.5 

4.5 

4.5 

63 

.8 

7.8 

4.7 

3.7 

69 

.8 

6.2 

5.37 

7.1 

81 

2 

20 

5.95 

6.7 

86 

1.3 

12 

5.55 

5.0 
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